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The optical properties of layered tungsten ditelluride have been measured over a wide temper¬ 
ature and frequency range for light polarized in the a-b planes. A striking low-frequency plasma 
edge develops in the reflectance at low temperature where this material is a perfectly compensated 
semimetal. The optical conductivity is described using a two-Drude model which treats the elec¬ 
tron and hole pockets as separate electronic subsystems. At low temperature, one scattering rate 
collapses by over two orders of magnitude, while the other also undergoes a significant, but less 
dramatic, decrease; both scattering rates appear to display the quadratic temperature dependence 
expected for a Fermi liquid. First principles electronic structure calculations reveal that the low- 
lying optical excitations are due to direct transitions between the bands associated with the electron 
and hole pockets. 

PACS numbers: 78.30.-j, 72.15.Lh, 71.20.Be 


Materials containing tellurium display a wide range of 
electronic properties that make them of interest for a va¬ 
riety of different applications, including photovoltaics [T] , 
thermoelectrics [2], p-type semiconductors [3], topologi¬ 
cal insulators |4j , and superconductors [5] . To add to this 
wide variety of phenomena it has recently been proposed 
that the layered transition metal dichalcogenide WTe 2 is 
a perfectly compensated semimetal at low temperature 
[5], i.e., the number of electrons and holes per unit vol¬ 
ume is identical. A simple two-band model predicts that 
the magnetoresistance will saturate at high fields, unless 
the material is compensated, in which case the magne¬ 
toresistance will increase quadratically with held [7] . The 
magnetoresistance of WTe 2 increases in precisely such a 
fashion; moreover, it becomes extremely large at low tem¬ 
perature and and shows no signs of saturation, even at 
very high helds [6]. Electronic structure and transport 
studies indicate that there are electron and hole pockets 
in this material mis], and angle resolved photoemission 
studies show that at low temperature they are almost ex¬ 
actly the same size [TOl [11] , resulting in a perfectly com¬ 
pensated semimetal. The zero-held behavior is equally 
dramatic with a very large residual resistivity ratio of 
p(300K)/p(5K) ~ 300, with p(5K) ~ 2 x 10"® flcm 
[5|. It is expected that such a dramatic change in the 
resistivity should produce striking changes in the optical 
properties of this material. 

In this study the optical properties of WTe 2 are deter¬ 
mined in the a-b plane in the absence of a magnetic held 
over a wide frequency and temperature range. The re- 
hectance reveals the formation of a striking low-frequency 
plasma edge at low temperature. The complex con¬ 
ductivity indicates that the electron and hole pockets 
give rise to carriers with dramatically different scatter¬ 
ing rates. While both follow a quadratic temperature 
dependence and may be thought of as Fermi liquids, one 
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Figure 1. (Color online) The temperature dependence of the 
reflectance of in the infrared region of WTe 2 for light polarized 
in the a-b planes, revealing the formation of a sharp plasma 
edge in the reflectance with decreasing temperature. Inset: 
The unit cell is shown for the b-c face projected along the a 
axis. 

scattering rate decreases by more than two orders of mag¬ 
nitude at low temperature and is responsible for the for¬ 
mation of a plasma edge in the reflectance. Previous 
electronic structure calculations are reproduced, showing 
electron and hole pockets of roughly equal size, and ex¬ 
tended to include the real part of the in-plane optical 
conductivity. The low-energy interband transitions are 
shown to originate from direct transitions between the 
bands associated with the hole and electron pockets. 

Single crystals of WTe 2 were grown using a bromine 
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Figure 2. (Color online) The temperature dependence of the 
real part of the optical conductivity of WTe 2 for light polar¬ 
ized in the a-b planes. The symbols along the y axis denote 
the values for aoc obtained from transport measurements [6] . 
Inset: The optical conductivity for a restricted set of temper¬ 
atures over a wider frequency range. 


vapor transport method [^. WTe 2 crystallizes in the 
orthorhombic Pmn2i space group |12] . where the tung¬ 
sten atoms form chains along the a axis and lie between 
sheets of tellurium atoms, forming the a-b planes. The 
sheets stack along the c axis and are easily exfoliated 
(inset of Fig. [^. The foil-like platelets were glued to a 
metal backing plate for support and then mounted on 
optically-black cones. The temperature dependence of 
the reflectance, shown in Fig. was measured for light 
polarized in the a-b planes over a wide frequency range 
(~ 2 meV to 3 eV) using an in situ overcoating technique 
[13j . The low-frequency reflectance is characteristic of a 
metallic response, increasing dramatically with decreas¬ 
ing temperature and revealing a sharp plasma edge at 
^ 300 — 400 cm“^. The optical conductivity has been 
determined from a Kramers-Kronig analysis of the re¬ 
flectance [14] , which requires that the reflectance be sup¬ 
plied over the entire frequency range; below the low¬ 
est measured frequency a metallic Hagen-Rubens form 
has been employed, 1 — i?(a;) oc ^/uj, while above the 
highest-measured frequency the reflectance is assumed 
to decrease as l/tu up to 1.5 x 10® cm“^, above which a 
free-electron l/w^ extrapolation was used [15]. 

The optical conductivity in the infrared region is shown 
in Fig. At high temperatures the response may be 
reasonably described by a single Drude component (a 
Lorentzian centered at zero frequency where the full 
width at half maximum is the free-carrier scattering 


rate), that gives way to a mid-infrared response that is 
dominated by interband transitions. As the temperature 
is reduced this Drude-like feature narrows considerably 
leading to a transfer of spectral weight (area under the 
conductivity curve) from high to low frequency, reveal¬ 
ing several interband transitions, shown more clearly in 
the inset of Fig. Above 2000 cm“^ (0.25 eV) the op¬ 
tical conductivity increases linearly with frequency until 
roughly 1 eV; we argue that this linear behavior is due 
to the superposition of several different interband tran¬ 
sitions and not due to a strong renormalization of the 
free-carrier scattering rate [Hin]. At high temperature 
there is reasonable agreement with (Ti(w — > 0) and the 
transport values for udc = 1/p- However, at low tem¬ 
peratures a single Drude component is not capable of 
describing both the extremely high values for ctdc and 
the shape of the low-frequency conductivity. Given that 
there are both electron and hole pockets in this mate¬ 
rial, it is appropriate to employ a model that accounts 
for both types of carriers [TH]. The complex dielectric 
function e = ei -f ie 2 of the so-called two-Drude model is 
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where Coo is the real part at high frequency. In the 
first sum ujpj = Anuje^/m* and l/r^ are the square of 
the plasma frequency and scattering rate for the delo¬ 
calized (Drude) carriers in the jth band, respectively, 
and Uj and m* are the carrier concentration and ef¬ 
fective mass. In the second summation, a;*,, 7 ^ and 
Dfc are the position, width, and strength of the fcth vi¬ 
bration or bound excitation. The complex conductiv¬ 
ity is cr(cLi) = (Ti -f ia 2 = —2TTiuj[e(u!) — eooj/^^o, where 
Zq — 377 n is the impedance of free space, yielding units 
for the conductivity of D“^cm“^. 

The two-Drude model is fit simultaneously to both cti 
and (72 using a non-linear least-squares method. The fit 
to the data at 100 K shown in Fig. [^indicates that the op¬ 
tical conductivity can be reproduced quite well using two 
Drude components, in addition to a series of Lorentzian 
oscillators at ~ 800, 1100 and 1950 cm“^ (~ 100, 136 and 
240 meV); other bound excitations above ~ 3500 cm“^ 
(430 meV) are also included in the total fit but not 
discussed. At 100 K the two Drude contributions are 
of roughly equal weight with ujp^i = 3000 ± 200 cm“^, 
1/ri = 26 ± 5 cm“^ and u!p ^2 = 3400 ± 200 cm“^, 
1 /t 2 = 89 ± 12 cm“^; fits at 80 and 150 K returned 
similar values for the plasma frequencies, so they have 
been assumed to be temperature independent and are 
not fitted parameters. While fitting to both scattering 
rates at 100 K works reasonably well, at higher tem¬ 
peratures the broad nature of the conductivity leads to 
larger values for the scattering rates and increasing un¬ 
certainties. At lower temperatures the dramatic increase 
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Figure 3. (Color online) The fit to the complex optical con¬ 
ductivity of WTe 2 at 100 K. The fit is composed of two 
Drude responses (dashed lines) and several Lorentzian fea¬ 
tures (dot-dash lines); the linear combination (long-dash line) 
reproduces the real part of the optical conductivity quite well. 
Inset: The fitted values for l/n (circles) and l/r 2 (squares); 
the dashed lines are drawn as a guide to the eye. 

in ctdc below 100 K suggests that 1/ti is becoming too 
small to be reliably determined from fits to d'{uj). In 
order to provide an anchor for the hts to the scatter¬ 
ing rates, 1/ri is first estimated from the Drude formula 
1/ri = 2TT ujp i/(Z q ao). For those temperatures where a 
direct comparison between transport and optics is pos¬ 
sible, it appears that (Tq = cri(a; —0) ~ 1.2 (Tdc, so 
we have used the approximation uo — 1-2 (Tdc- Treating 
1 /ti as static, only 1 is fit to the optical conductivity, 
allowing the contribution cro 2 = 27r Wp_ 2 '’’ 2/^0 to be cal¬ 
culated where ctq = coi + cro 2 ; the iterative improvement 
in 1/ri is determined from 1/ri = 2t: lUp i/{Z q am)- The 
temperature dependence of the transport values of p, ag, 
and the scattering rates determined in this fashion are 
summarized in Table |T] The temperature dependence of 
1/ri and 1 /t 2 is shown in the inset of Fig. as circles and 
squares, respectively. While there is some uncertainty of 
the scattering rates at high temperature, in both cases 
they follow a roughly quadratic temperature dependence, 
suggesting both sets of carriers resemble a Fermi liquid 
[IMB- While l/r 2 decreases by roughly an order of 
magnitude at low temperature, it does not rival the dra¬ 
matic collapse of 1/ri from ~ 110 cm“^ at 295 K to 
~ 0.3 cm“^ at 5 K, associated with the formation of the 
plasma edge in the reflectance. At low temperature the 
extremely large, non-saturating magnetoresistance indi¬ 
cates that this material is perfectly compensated m, 


Table I. Transport values p and ctq, and optical scattering 
rates 1/ri and l/r 2 . 


T 

(K) 

P’" 

{pQ cm) 

ao 

(x 10^ D“^cm“^) 

l/ri 

(cm-i) 

1 / 7-2 

(cm-i) 

295 

580 

2.07 

110 ±22 

320 ± 45 

200 

400 

3.00 

84± 16 

160 ± 22 

150 

280 

4.28 

60± 12 

109 ± 15 

100 

150 

8.00 

26 ±5 

89± 12 

80 

100 

12.0 

16 ±3 

80± 11 

60 

70 

17.2 

12 ±2 

47 ±6 

40 

30 

40.0 

4.1 ± 0.8 

54 ±7 

20 

7 

171 

0.90 ±0.2 

38 ±5 

6 

2 

600 

0.25 ±0.05 

34 ±4 


Ref. |6] 


i.e. ni = 77,2. From the values of LUpp and ujp^ 2 , it may 
be shown that carrier masses are similar, m 2 — O.Sttii, in 
good agreement with Shubnikov-de-Hass oscillation mea¬ 
surements [22] • 

There are a series of excitations in the infrared at 
~ 100,136 and 240 meV, as well as several others above 
0.4 eV; these absorptions are in general quite broad and 
have a substantial amount of spectral weight at low fre¬ 
quency. Previous band structure studies of WTe 2 re¬ 
vealed a semimetallic character with small electron and 
hole pockets 00; these results have been reproduced 
and extended in this work (the details of the calculations 
[ 2 M 25 ] are discussed in the supplementary materials). 

The calculated electronic band structure including the 
effects of spin-orbit coupling is shown in Fig. Ba) for 
paths between several different high-symmetry points in 
the orthorhombic Brillouin zone, revealing both electron 
and hole bands in the F — X direction, in good agreement 
with previous results 00 . A more detailed examination 
of the band structure, (supplementary material), reveals 
several bands that lie just above and below the Fermi 
level ep that may be responsible for the low-energy in¬ 
terband transitions. To explore this further the total 
density of states (DOS) for the W and Te atoms is cal¬ 
culated using a fine fc-point mesh (10 000 k points) and 
shown in Fig. Bh)- The low value for the DOS at ep is 
consistent with the semimetallic nature of this material. 
The prominent increase in the DOS in the 200 — 500 meV 
above ep is larger than 50 — 150 meV interval over which 
interband transitions are observed in the optical conduc¬ 
tivity. However, the optical conductivity is a reflection 
of the joint density of states (JDOS) rather than just 
the DOS. The real part of the optical conductivity has 
been calculated |26j using the same fine /c-point mesh 
for a^^x and ay^y (a and b axes, respectively), shown in 
Fig. B c), and is compared to the experimental value for 
the conductivity at 5 K with the Drude contributions 
removed. Both ax^x and ay^y have a very weak absorp¬ 
tion at about 25 meV with the onset of much stronger 
absorptions occurring above about 50 meV, with ax,x dis- 
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Figure 4. (Color online) (a) The calculated electronic struc¬ 
ture WTe 2 including the effects of spin-orbit coupling between 
several different high symmetry points of the orthorhombic 
Brillouin zone, (b) The total density of states (DOS), includ¬ 
ing the contributions from the W and Te atoms, (c) The real 
part of the calculated optical conductivity ay^y and aver¬ 
age of the two, compared to the measured optical conductivity 
at 5 K with the Drude terms removed. 

playing prominent features at ~ 80 and 160 meV, while 
only the 160 meV peak is present cry,y', above ~ 200 meV 
the two conductivities continue to increase steadily with 
frequency. The average of the two curves reproduces the 
experimental data reasonably well. The details of the 
contributions to the optical conductivity, (supplementary 
material), reveal that these two strong absorptions at 80 
and 160 meV originate from direct transitions between 
the bands associated with the hole and electron pock¬ 
ets. The feature at 60 meV, which occurs at too high 
an energy to be a lattice vibration m , remains unex¬ 
plained; however, it has been shown that there are very 
low-energy excitations in this material. Because this fea¬ 
ture only emerges at low temperature where the material 
is perfectly compensated, it is likely that it is due to a 
low-energy interband transition that is sensitive upon the 
details of the band structure that may not captured in 
this calculation. 

To conclude, the complex optical properties of the per¬ 
fectly compensated semimetal WTe 2 have been measured 
in the a-b planes for a variety of temperatures over a wide 


frequency range. The electron and hole pockets are de¬ 
scribed using the two-Drude model, revealing that as the 
temperature is reduced one scattering rate collapses by 
over two orders of magnitude, resulting in the formation 
of a striking plasma edge in the reflectance, while the 
other scattering rate is considerably larger and does not 
decrease by nearly the same amount. Both sets of carriers 
may be described as Fermi liquids. The calculated op¬ 
tical conductivity indicates that low-energy absorptions 
are due to direct transitions between the bands associ¬ 
ated with the hole and electron pockets. The emergence 
of a narrow low-energy absorption at low temperature is 
thought to be a low-lying interband transition related to 
subtle changes in the electronic structure that result in a 
perfectly compensated semimetal. 
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ELECTRONIC STRUCTURE CALCULATIONS 


Tungsten ditelluride, WTe 2 , crystallizes in the orthorhombic Pmn2i space group [1], where the tungsten atoms 
form chains along the a axis and lie between sheets of tellurium atoms, forming the a-b planes; the unit cell 
is reproduced from the main text in Fig. SI for convenience. The electronic structure of WTe 2 was calculated 
using both the local spin density approximation (LSDA) as well as the generalized gra¬ 
dient approximation (GGA) using the full-potential linearized augmented plane-wave 
(FP-LAPW) method [2] with local-orbital extensions [3] in the WIEN2k implemen¬ 
tation [4]. The system was assumed to be non-magnetic. An examination of different 
Monkhorst-Pack fc-point meshes indicated that a 9 x 5 x 2 mesh and Rmtkmax = 7.75 
was sufficient for good energy convergence. The geometry of the unit cell was refined 
through an iterative process whereby the volume was optimized with respect to the 
total energy while the c/a ratio remained fixed. The atomic fractional coordinates 
were then relaxed with respect to the total force, typically resulting in residual forces 
of less than 0.2 mRy/a.u. For both of these procedures spin-orbit coupling is ignored. 

This approach was repeated until no further improvement was obtained. A compar¬ 
ison of the experimental and calculated (relaxed) unit cell parameters are shown in 
Table. I. 

The electronic band structure has been calculated with LSDA and spin-orbit cou¬ 
pling for several different paths between high-symmetry points in the orthorhombic 
unit cell, shown in Fig. S2. This calculation accurately reproduces the work of Ali et 
al. [5], as well as an earlier work [6]. The relatively small number of bands crossing 
the Fermi level (ep) illustrates the semimetallic character of WTe 2 . In particular, 
there is an electron-like band and a hole-like band crossing Fermi surface along the 
F — X direction. Figure S2 also reveals that there are several other points in the 
Brillouin zone at which hole and electron bands just manage to cross the Fermi sur¬ 
face; however, the electronic structure is quite sensitive to the unit cell geometry, and 
that even a small increase (or decrease) in the Fermi level might move these bands 
above or below ep. While the paths along the high-symmetry directions are useful, 
the contributions of the hole and electron bands to the transport is best understood 



Figure SI. (Color online) The 
unit cell is shown for the b-c 
face projected along the a axis. 
The tungsten atoms are sur¬ 
rounded by tellurium octahe- 
dra that are distorted due to 
the zig-zag nature of the tung¬ 
sten chains. 


Table I. The experimental and theoretical lattice constants and atomic fractional coordinates for the relaxed structure of WTe 2 
in the orthorhombic Pmn2\ space group. 


a(A) 

Experiment 

a 

Theory (LSDA) 

Theory (GGA) 


3.5306 



3.4657 



3.5423 


6(A) 


6.3443 



6.2274 



6.3653 


cU) 

W1 {\xy) 

0.39890 

14.2123 

0.50030 

0.39866 

13.9509 

0.50076 

0.39859 

14.2595 

0.50039 

W2 (Oxy) 

0.04310 


0.48490 

0.04368 


0.48491 

0.04393 


0.48502 

Tel i^xy) 

0.79410 


0.59540 

0.79507 


0.59548 

0.79470 


0.59428 

Tel (Oxy) 

0.30040 


0.63990 

0.29920 


0.64169 

0.30052 


0.63948 

Tel {Oxy) 

0.85620 


0.84510 

0.85526 


0.84396 

0.85649 


0.84588 

Tel i^xy) 

0.35180 


0.89020 

0.35271 


0.89030 

0.35225 


0.89129 


Ref. 1. 
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Figure S2. (Color online) The calculated LSDA electronic band structure of WTe 2 including the effects of spin-orbit coupling 
shown for several different paths between high-symmetry points for an orthorhombic Brillouin zone (generated using 500 k 
points). 


from the nature of the Fermi surface. 


(a) 


- - 

(b) 

H 

(c) 

it ig 

(d) 

• * • • 



Figure S3. (Color online) The (a) small and (b) large hole LSDA Fermi surfaces, (c) Large and (d) small electron Fermi 
surfaces, (e) The merged Fermi surface of the electron (yellow) and hole (blue) pockets in the orthorhombic Brillouin zone 
along the F — X direction. 
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FERMI SURFACE 


The Fermi surface has been calculated using LSDA and a dense fc-point mesh (10 000 k points, resulting in a 
41 X 23 X 10 mesh) [7]. The resulting Fermi surface is shown in Fig. S3 and consists of a total of four pockets; two 
hole pockets, one small and one large, shown in Figs. S3(a) and (b), respectively, and two electron pockets, one large 
and one small, shown in Figs. S3(c) and (d), respectively. The merged Fermi surface shown in Fig. S3(e) consists of 
the electron (yellow) and hole (blue) pockets slightly displaced from the origin of the Brillouin zone along the F — X 
direction, in agreement with a previous calculation [5] and several experimental observations [8, 9]. 


DIELECTRIC TENSOR 


In the main text the real part of the optical conductivity has been calculated [10] from the the imaginary part of 
the dielectric function, ax,x = 27ra; 5 ex,x/Zo, using a fine k point mesh (10 000 k points, yielding a 41 x 23 x 10 mesh); 
Zq ~ 377 fl is the impedance of free space, resulting the units for the conductivity of n“^cm“^. 

In this case, the free-carrier contribution is not calculated, so that the imaginary part of the dielectric function is 
the sum of the contributions to the dielectric tensor over all the allowed interband transitions. There are a total of 
617 energy bands, with the last occupied band occurring at 276. Only direct (q = 0) transitions are considered from 
occupied to unoccupied states; this being the case, the bands near the Fermi surface will be of particular importance as 
they will make the most significant contributions to the imaginary part of the dielectric function and thus the optical 


conductivity. The energy bands within 1 eV of the Fermi 
Note that in the main text both and ay^y have 

been calculated; however, whereas ax,x displays promi¬ 
nent structures at 650 cm“^ (80 meV) and 1300 cm“^ 
(160 meV), the first prominent feature in ay^y only oc¬ 
curs at about 1300 cm“^, therefore this calculation will 
focus on the result for imaginary part of [H]- Tran¬ 
sitions between all energy bands have calculated; how¬ 
ever, we will only consider contributions to 3ex,x ^ 0.5. 
These transitions are shown in Fig. S5. For bands 235 
and 236, there are only weak contributions from tran¬ 
sitions to bands 241 and 242 that only occur at ener¬ 
gies above about 3400 cm“^ (420 meV); however, for 
bands 237 and 238, we see contributions from an in¬ 
creasing number of bands (239 —>■ 244) with a much 
lower onset of about 400 cm“^ (50 meV) with noticeable 
peaks at about 700 (87 meV) and 2500 cm“^ (310 meV). 
The contribution from band 239 is particulary interest¬ 
ing as there is an unusually low-energy transition at about 
150 cm“^ (19 meV) with another prominent absorption at 
1400 cm“^ (174 meV); for band 240 the strongest absorp¬ 
tion is at about 800 cm“^ (99 meV). For bands 241 and 
242, the remaining contributions are fairly small. Thus, 
the two strongest contributions are from the transitions 
lated to the electron and hole pockets that form the Fern 
transitions originate from these energy bands. 


level along the F — X direction have been labeled in Fig S4. 



Figure S4. (Color online) The indexed LSDA energy bands 
near the Fermi level along the F — X direction in the or¬ 
thorhombic Brillouin zone. 

239 —>■ 242 and 240 —)■ 241, which are the energy bands re- 
h surface so we may conclude that the low-energy interband 


The sum of all of the contributions to the imaginary part of the complex dielectric function 3 ex,x is shown in 
Fig. S6, with the calculated real part of the optical conductivity shown in the inset [note this is same quantity that 
is shown in Fig. 4(c)]. The peak at about 150 cm“^, while prominent in the imaginary part of Cx^x, is diminished in 
the optical conductivity due to the effects of frequency scaling. 
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Figure S5. (Color online) The LSDA contributions to the imaginary part of the dielectric function Syea;,^; for the various 
transitions between occupied {bi) to unoccupied (6/) bands, beginning at 235 and ending at 242. In each panel, the initial 
band is above the horizontal line, while the final band (or bands) is listed below the line; note that the choice of bands was 
determined from the criteria at the contribution should be larger than ~ 0.5. 
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Figure S6. (Color online) The imaginary part of the dielectric function from a summation of all of the individual interband 
contributions, determined using LSDA (solid line) and GGA (dashed line) with spin-orbit coupling. Inset: The real part of the 
optical conductivity ax,x. calculated from the imaginary part of complex dielectric function ex,x using LSDA (solid line) and 
GGA (dashed line). 
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